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^SJ , Abstract. Age constraints are most often placed on globular clusters by comparing their CMDs 
with theoretical isochrones. The recent discoveries of detached, eclipsing binaries in such systems 
by the Cluster AgeS Experiment (CASE) provide new insights into their ages and, at the same 
time, much needed tests of stellar evolution models. We describe efforts to model the properties 
of the detached, eclipsing binary V69 in 47 Tuc and compare age constraints derived from stellar 
(-H ' evolution models of V69A and B with ages obtained from fitting isochrones to the cluster CMD. 

We determine, under reasonable assumptions of distance, reddening, and metallicity, the extent 
I ' to which it is possible to simultaneously constrain the age and He content of 47 Tuc. 
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J> 1. Background and Previous Results 

' 47 Tuc is among the closest and most carefully studied globular clusters in the Galaxy. Table 

' [T] lists the best available estimates of [Fe/H], [a/Fe], distance, and reddening along with a 

, representative sample of ages from recent studies. 

\6 ' 

^ • 2. Stellar Evolution Models 

. The stellar evolution models utilized in this contribution were computed using the Dartmouth 

' Stellar Evolution Program (DSEP; Dotter et al. 2007, 2008). The models include the effects of 

, partially inhibited microscopic diffusion of He and metals (Chaboyer et al. 2001). The models 
were converted to the observational plane using PHOENIX synthetic spectra. Stellar evolution 

$H tracks and isochrones were computed for [Fe/H]= -0.8, -0.75, and -0.7; [a/Fe]= 0, -1-0.2, and 

p^. , +0.4; and Y=0.24, 0.255, 0.27, 0.285, and 0.3. The color transformations include the effects of 
a-enhancement but all assume Y~0.25. 



3. Age Constraints from the CMD 

The analysis performed in this section uses the color-magnitude diagram (CMD) of 47 Tuc 
from the ACS Survey of Galactic Globular Clusters (Sarajedini et al. 2007; Anderson et al. 
2008). The data have been culled (by removing stars with large photometric errors, see §7 of 
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Table 1. 47 Tuc: Basic Parameters and Previous Age Results 



Parameter 


Value 


Source 


[Fc/H] 


-0.76±0.01±0.04 
-0.75±0.01±0.04 


Koch & McWilliam (2008) 
McWilliam & Bernstein (2008) 


[a/Fe] 


~0.4 
~0.3 


Koch & McWilliam (2008) 
McWilliam & Bernstein (2008) 


DMv 


13.35±0.08 


Thompson et al. (2009) 


E(B-V) 


0.0320±0.0004 


Schlegel et al. (1998) 


Age (Gyr) 


10-13 
~11.3 
11.0±1.4 
10.7±1.0 
11.5±0.8 
12.5±1.5 


Salaris et al. (2007) 
Gratton et al. (2003) 
Percival et al. (2002) 
Salaris & Weiss (2002) 

VandenBerg (2000) 
Liu & Chaboyer (2000) 



Table 2. Parameters of V69 from Thompson et al. (2009) 



1 Parameter 


Primary 


Secondary | 


1 M/Mq 


0.8762±0.0048 


0.8588±0.0060 | 


1 R/ R@ 


1.3148±0.0051 


1.1616±0.0062 1 


1 L/Lq 


1.94±0.21 


1.53±0.17 1 



Anderson et al. 2008) in order to more clearly delineate the main sequence but the final CMD 
still contains ~50,000 stars. 

In order to simplify the analysis, the apparent distance modulus and reddening were fixed to 
those values listed in Table [1] Hence the only uncertainties in the ages derived from isochrone 
fitting are due to the inherent scatter in the data and the allowed range of metallicities. 

The basic results from fitting isochrones to the CMD are as follows: (i) the main sequence and 
red giant branch impose little or no constraint on Y variations at the ~0.05 level (see Figure[T|; 
(ii) the models favor a level of [a/Fe] that is lower than the spectroscopic value (see Figure [2| 
and indeed appear to rule out [a/Fe] =+0.4; and (iii) the range of possible [Fe/H] values listed 
in Table 1 all give acceptable fits to the CMD assuming the preferred value of [a/Fe] is used 
(see Figure [3}. 

Figure |4] shows how age and [Fe/H] (assuming [a/Fe]=+0.2 and Y~0.25) are correlated. The 
solid line is the best fit value as a function of [Fe/H] and the dashed lines represent the error bars 
that arise only from the finite width of the subgiant branch. A value of Y=0.255 was adopted for 
the plot but any Y value within the assumed range would give a similar result. The isochrone fits 
to the CMD yield an age of 11.5±0.75 Gyr where the uncertainty due to the fitting procedure 
is ~0.3 Gyr and the rest is due to ~0.05 dex uncertainty in [Fe/H]. 

4. Age Constraints from the Binary V69 

A careful analysis by Thompson et al. (2009) of spectroscopic radial velocity measurements 
and photometric light curves provides high precision estimates of the fundamental parameters of 
V69, see Table [2] While the masses and radii are measured to better than 1%, the luminosities 
are measured to no better than ~10% at present (though additional near- infrared light curves 
will substantially improve this number). 

In order to measure the age of 47 Tuc using V69, individual stellar evolution tracks were 
computed for six diff'erent masses: three for each star, encompassing the central value and quoted 
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uncertainties in Table (2] These models were calculated for the same composition grid as the 
isochrones described in section 2. The mass-radius (M-R) and mass-luminosity (M-L) relations 
in the tracks were used to determine the range of age for which each track lay within the error 
bars for the corresponding parameter (R or L) of each component. We make two important 
general comments: (i) the ages derived from each relation are in excellent agreement between 
the two stars and (ii) the agreement between the two methods in a given star are quite sensitive 
to Y. For either star, the maximum Y value for which the M-R and M-L results overlap increases 
with [Fe/H]. 

The results of the V69 age analysis are presented in the Figure [5] The solid lines indicate how 
the measured age varies with [Fe/H] for the Y value listed on the figure. The dashed lines above 
and below the line for Y=0.255 indicate the combined uncertainties in M, R, and L, with ^0.5 
Gyr due to the luminosity uncertainty alone. Figure [5] clearly shows that the age-[Fe/H] relation 
is complimentary to that of the CMD method (Figure U)). Of equal importance is the fact that 
since the masses of the stars are known, there is a significant sensitivity to the He abundance 
that is lacking in the CMD method. 

5. Combined Constraints 

It was demonstrated in §3 that the CMD age- [Fe/H] relation (assuming that both the distance 
and reddening are well-constrained) is linear for -0.8 < [Fc/H] < -0.7 and that age decreases 
as [Fe/H] increases. Conversely, the result from §4 showed that while the age- [Fe/H] relation 
derived from V69 was also linear over the range of [Fe/H] considered, the age increased as 
[Fe/H] increased. Combining the two methods is therefore quite powerful because the slopes 
have opposite signs. The region where the two overlap has a finite size and gives preferred 
values for age, [Fe/H], and Y simultaneously, though ultimately the size of the overlap region is 
dominated by the uncertainties in both methods. 

Figure [S] presents the results given in Figures U and [S] together so that the situation described 
in the preceding paragraph can be clearly seen. The figure indicates that the binary and CMD 
methods agree best if the He abundance is low (Y 0.24) while a value of Y 0.25 would 
have been more preferable given the constraints from WMAP and Big Bang nucleosynthesis, 
see Spergel et al. (2003). 

If the bias towards [a/Fe] = -1-0.2 derived from isochrone fits to the CMD in section 2 is 
relaxed and the spectroscopic value (-1-0.4) is adopted instead, then the combined results favor 
Y ~ 0.25. However, it is difficult to reconcile the larger [a/Fe] ratio given the precise reddening. 

Further improvements to the distance and luminosity of V69 from near infrared light curves 
will significantly reduce the uncertainties and provide a much more stringent test of the stellar 
evolution models' ability to simultaneously fit the components of V69 and yield an age consistent 
with the turnoff age. As such, the ongoing CASE project represents the most exacting test of 
stellar evolution models in globular clusters to date and is certain to provide invaluable insights 
into the ages-and other aspects-of such systems. 
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Figure 1. The effect of varying Y at fixed [Fe/H] and [a/Fe]. Shown are isochrones with 

Y=0.24, 0.255, 0.27, and 0.285. 
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Figure 2. The effect of varying [a/Fe] at fixed [Fe/H] and Y. 
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Figure 3. The effect of varying [Fe/H] at fixed [a/Fe] and Y. 
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Figure 4. Isochrone fitting results in tiio age [Fe/H] plane. The solid line represents the best fits 
while the dashed lines indicate the uncertainties due only to the width of the subgiant branch. 
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Figure 5. V69 results in tiie age-[Fe/H] plane. The solid lines correspond to results from the 
labeled Y values. The dashed lines about the Y=0.255 line indicate the errors for that particular 
case which are typical of all cases. 
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Figure 6. Combined isochrone-CMD and V69 results in the age-[Fe/H] plane. 



